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Modified quark fragmentation functions in deeply inelastic eA collisions and their QCD evolution 
equations are derived for the first time in the framework of multiple parton scattering. Induced 
radiation gives rise to additional terms in the evolution equations and thus softens the modified 
quark fragmentation functions. The results in the next-leading-twist depend on both diagonal and 
off-diagonal twist-four parton distributions and the combination of which clearly manifests the LPM 
interference pattern. The predicted modification depends quadratically on the nuclear size {A'^^'^). 
Generalization to the case of hot QCD medium is also discussed. 
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The propagation of an energetic parton and its induced 
energy loss has been proposed as a probe of the proper- 
ties of dense matter formed in high-energy nuclear colli- 
sions [0,11 . Recent theoretical studies |3|-||] show that a 
fast parton will lose a significant amount of energy via 
induced radiation when it propagates through a hot par- 
tonic matter. One cannot directly measure the energy 
loss of partons because they are not final experimentally 
observed particles. However, parton energy loss does 
lead to modification of the final particle spectra. There- 
fore, one can only study the parton energy loss indirectly 
by measuring the modification of the parton fragmenta- 
tion functions in semi-inclusive processes like eA or 7-jet 
events in AA collisions [7\ or the inclusive spectra at large 
transverse momentum 1 2 ,P . 

In this Letter, we report our first study and derivation 
of the QCD evolution equations for the medium-modified 
fragmentation functions in the simplest case of deeply in- 
elastic eA scattering (DIS). The induced gluon radiation 
due to multiple parton scattering gives rise to additional 
terms in the modified QCD evolution equations that 
soften the modified fragmentation functions. Utilizing 
the generalized factorization of higher-twist (HT) parton 
distributions ||^ , we show that these additional HT terms 
depend on both the diagonal and off-diagonal twist-four 
parton distributions, the combination of which clearly 
manifests the Landau-Migdal-Pomeranchuk (LPM) in- 
terference pattern. Using estimates of these twist-four 
parton matrix elements from other processes such as the 
Pt broadening of Drell-Yan dilepton in pA collisions, we 
predict the modification of the effective quark fragmen- 
tation functions and their dependence on the parton en- 
ergy and nuclear size. We also estimate the quark energy 
loss defined as the total energy carried by gluons from 
induced radiation. 

We consider the following semi-inclusive process in 
the deeply inelastic eA scattering, e(Li) -I- A{p) — > 
e{L2) + h{£h)+X , where Li and L2 are the four momenta 
of the incoming and the outgoing leptons, £h is the ob- 



served hadron momentum, p and q = L2 — Li denoted as 
p — [p^ ,0,0±], q = [—Q'^/2q^,q^,0±], are the momen- 
tum per nucleon in the nucleus with the atomic number 
A and the momentum transfer, respectively. The differ- 
ential cross section for the semi-inclusive process can be 
expressed as 
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where s = (p + Li)'^ and oem is the electromagnetic 
(EM) coupling constant. The leptonic tensor is given by 
L^i, = l/2Tr(7 • _Li7^7 • i27i/) while the semi-inclusive 
hadronic tensor is defined as, 
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X 

X 2Tr5'^{q + p~px -4) 



(2) 



where runs over all possible final states and — 
J2q Gq^'qltii^q IS the hadrouic EM current. 

In the parton model with coUinear factorization ap- 
proximation and to the leading-twist (LT) the semi- 
inclusive cross section factorizes into a product of par- 
ton distributions, parton fragmentation functions and the 
partonic cross section. Therefore, to the leading order in 
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where the momentum fraction carried by the hadron is 
defined as Zh = ^h/'l~ ^^id xb = Q^/2p^q~ is the 
Bjorken variable. /Zj and /i^ are the factorization scales 
for the initial quark distributions f^[x,y?j) in a nucleus 
and the fragmentation functions Dq^h{zh, n^), respec- 
tively. Including all leading log radiative corrections, the 
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renormalized quark fragmentation function Dq^h{zhi A^^) 
satisfies the QCD evolution equation ||ic| ]. 
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FIG. 1. Diagrams for rescattering with gluons (a) and 
quarks (b) without and with gluon radiation (c) in deeply in- 
elastic eA scattering. Possible cuts are shown by the dashed 
lines 

In this paper we will consider contributions of quark 
rescattering with partons from another nucleon inside the 
nucleus. Such contributions are proportional to the nu- 
clear size A^/'^ For large enough A, we can neglect 
other A-independent HT effects. For large in DIS, 
it may suffice to only consider one rescattering. The 
contributions of one rescattering can be treated as HT 
corrections to the LT results. We work in a frame- 
work p2| in which the twist-four contributions can be 
expressed as the convolution of the partonic hard parts 
and four-parton matrix elements j^]. At the lowest order, 
rescattering without gluon radiation as shown in Fig. [l|(a) 
broadens the transverse momentum of the leading jet ]11|] 
but contribute little to parton energy loss. One can also 
neglect rescattering with another quark in Fig. 0(b). 

The dominant HT contributions to the QCD evolu- 
tion of the fragmentation functions come from radiative 
processes involving rescattering with a gluon from an- 
other nucleon as illustrated by the central-cut diagram in 
Fig. [^(c). Kinematics only allows two poles, one at each 
side of the central-cut, out of the four propagators in the 
diagram. This leads to four possible combinations each 
give different momentum fractions to the initial partons. 
In one case, the initial gluon has X2 = xl + xd which is 
finite when kx 0, where 
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It is the transverse momentum of the radiated gluon, kx 
is the initial gluon's intrinsic transverse momentum, and 
z — tq jq^ is the momentum fraction carried by the final 
quark. This corresponds to gluon radiation induced by 
the rescattering and is referred to as a double-hard pro- 
cess. In another combination, X2 — xjj which vanishes 
when fcT — > 0. In this case the rescattering is soft and 
the gluon radiation is induced by the initial hard photon- 
quark scattering. Such a process is called hard-soft. The 
four contributions from Fig. |l|(c) correspond to these two 
distinct processes and their interferences. Their sum has 
the form, 



H^Ji) (X (1 - e~"^f y-^ )(i - e^^^p -y^ )) 
X e"''p^'^y^-y^\ (5) 



This clearly manifests the LPM interference pattern 
caused by the destructive interferences between hard-soft 
and double- hard processes. The interference pattern is 
dictated by the gluon's formation time, r/ = 1/xlP~^, rel- 
ative to the nuclear size. The two processes completely 
cancel each other in the collinear limit when £t —> 0- 
Diagrams involving three-gluon vertices have exactly the 
same structure as Fig. |l|(c) , except that they have differ- 
ent momentum dependence and color factor in the hard 
part. 

We have considered all together 23 possible cut dia- 
grams, 14 of them are interferences between no and dou- 
ble rescattering (shown as the left and right-cut diagrams 
in Fig. |l|(c)) which cancel some of the contributions from 
central-cut diagrams. Including virtual corrections, we 
obtain |l^] the leading HT contribution from rescatter- 
ing processes. 
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where 



Tq,{x,XL) = / — 
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is quark-gluon correlation function which essentially con- 
tains four independent four-parton matrix elements in a 
nucleus and xt — {k^)/2p+q^ = XB{k^)/Q^- With the 
definition of the -I- functions , the term proportional 
to the i5-function accounts for virtual corrections and 
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dz- 



{1^z^)T^^{x,xl)]. (8) 



One can similarly get the contribution from the gluon 
fragmentation. We will neglect the radiative corrections 
to processes such as Fig. ^b) that involve rescattering 
with a quark in the leading log approximation, because 
they can be shown to be proportional to |^ as 

compared to l/t^ in Eq. (^). 

Summing up all the leading contributions from LT 
and HT processes, we can effectively define the modified 
quark fragmentation function as 



dxfq(x, fij)Hlj^J{x,p, q)Dq^h(zh, 



(9) 

where for completeness should also include the 

HT contributions as studied by Mueller and Qiu [ p^ . 
The modified quark fragmentation function satisfies the 
following evolution equation 



To estimate the twist-four parton matrix elements, 
we generalize the approach by |^ to include the off- 
diagonal matrix elements. Assuming a Gaussian nuclear 
distribution in the rest frame, p{r) ^ exp(— r^/2i?^), 
Ra = 1.12A^/^ fm, we express in terms of single 
parton distributions, 

C 

T^ ix, xl) = — [fq{x){xT + xl)G{xt + xl) 

XA 

+ f^ix + xl)xtG{xt)] (1 ~ e-'^i^^"^), (14) 

where G{x) is the gluon distribution per nucleon in a 
nucleus, xa = 1/MRa, and M is the nucleon's mass. 
The off-diagonal terms involves transferring momentum 
Xl between different nucleons inside a nucleus and thus 
should be suppressed for large nuclear size or large mo- 
mentum fraction xl- Notice that Tf — 1/xlP'^ is the 
gluon's formation time. Thus, xl/xa — La/tj with 
La — RaM/p'^ being the nuclear size in our chosen 
frame. 
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with the modified splitting functions defined as 
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+ <5(l-z)Ar,^^(x,M')] , 
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where 'yq^qg{z) is the normal splitting functions [Q. We 
assume in the leading order that the gluon fragmentation 
function follows the normal QCD evolution equations. 

Solving the above equation is equivalent to summing 
all leading log twist-four contributions. As an approxi- 
mation, one can write the solution as, 

Dq^h{zh, ^i"^) = Dq^h{zh, fl'^) + ADq^h{zh, fJ."^) 



FIG. 2. The predicted modification to the quark fragmen- 
tation functions for three different values of initial quark en- 
ergy q~ = /2p'^XB- XA = 0.04 corresponds to A ~ 200 

Using the above approximation in Eq. ( |l^ ) and replac- 
ing the momentum fraction xl in the parton distributions 
by its average value (xl) ^ xa, we have 



ADq^h{zh,Q') 



2^ .-^ ^Aa'^ ^ 



Adq^hizh, XB , XA, Q"^), 



ADq^,{z,,t^) = ^ ^ - 
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^ dz 



X [A'Jq^qg{z, X, XLj^)Dq^h{zh/z, ^i^) 
-I- A'-fq^gq{z,X,XL,tT)Dg^h{zh/z,l-i'^)\ , 

(13) 

where Da^h{zh, are the normal fragmentation func- 
tions. Notice that there is no collinear divergence in the 
above integration because of the LPM effect in T^. Be- 
cause A7 is proportional to l/£^, ADq^h is suppressed 
by 1//^^ relative to the LT fragmentation function Dq^h- 



C = C[xtG{xt) + [xT + xa)G{xt + xa)], 

(15) 

where we choose the factorization scale /i^ = Q^. Shown 
in Fig. 2 are the numerical results of Ad{zh, xb,xa,Q^) 
for three different values of xb- The parameterization in 
Ref. of the normal fragmentation functions is used 
in our calculation. As we see from the numerical results, 
the modification to the shape of fragmentation function 
increases for larger values of corresponding to smaller 
quark energy q~ at fixed Q^. The modification also in- 
creases for smaller values of xa corresponding to larger 
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nuclear size. As shown in Eq. (y_5|), the magnitude of the 
modification depends quadratically on the nuclear size. 
This is because the LPM effect modifies the transverse 
momentum spectra of the radiated gluon such that the 
phase space for the momentum integral is limited. That 
limited phase space is proportional to the nuclear size. 
Together with the linear dependence of the parton cor- 
relation function, this leads to non-linear dependence of 
the total energy loss on the nuclear size. 

We can define the quark energy loss as the momentum 
carried by the emitted gluons. From Eq. (p^, we have, 

Jo *T .^0 ^'^ 



lO "'0 

CaCHs C Xb 



61n(-^) , {XA^XB < 1). (16) 
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According to Ref. p^ , p^ , C can be related to the trans- 
verse momentum broadening of Drell-Yan dilepton in pA 
collisions, 



A(/c|)dy 



27ras C 



0.022yl^/^ GeV^ 



(17) 



In the rest frame of the nucleus, the total energy loss is 
AE ^ qo{Azg) w 0.35 asA^/^ In GeV , (18) 

2XB 

which depends quadratically on the nuclear size. For 
XB = 0.1, as = 0.3 and A = 200, AE w 5.8 GeV which 
is consistent with the estimate in Ref [||. 

In summary, we have derived for the first time the 
evolution equations for the nuclear modified quark frag- 
mentation functions due to multiple scattering and in- 
duced radiation. We can generalize our results to the 
case of hot QCD matter. In our approximation, Eq. (14), 



the quark distribution which represents the initial quark 
production probability docs not enter into the modified 
evolution equations. What remains is the gluon matrix 
element. In the case of a hot QCD medium, the coor- 
dinates of the two gluon fields are not restricted to the 
nucleon size due to deconfinement in a nucleus. Rather, 
the final results will be sensitive to the correlation length 
of the hot medium jl9j . If there is a dramatic change of 
the correlation length during the QCD phase transition, 
one should also see a big change in the modification of 
the parton fragmentation functions. 



ACKNOWLEDGEMENTS 

We would like to thank J. Jalilian-Marian, J. Qiu 
and G. Sterman for helpful discussions. This work 
was supported by DOE under Contract No. DE-AC03- 
76SF00098 and DE-FG-02-96ER40989 and in part by 
NSFC under project 19928511. 



[2] 
[3] 

[4] 

[5] 

[6] 
[7] 

[8] 
[9] 

[10] 

[11] 
[12] 

[13] 
[14] 
[15] 
[16] 

[17] 
[18] 
[19] 



M. Gyulassy and M. Pliimer, Phys. Lett. B243, 432 
(1990); M. GyulassyM. Plumer, M. H. Thoma and X.-N. 
Wang, Nucl. Phys. A 538 (1992) 

X.-N. Wang and M. Gyulassy, Phys. Rev. Lett. 68, 1480 
(1992). 

M. Gyulassy and X.-N. Wang, Nucl. Phys. B420, 583 
(1994); X.-N. Wang, M. Gyulassy and M. Plumer, Phys. 
Rev. D 51, 3436 (1995). 

R. Baier, Yu. L. Dokshitzer, S. Peigne and D. Schiff, 
Phys. Lett. B345, 277 (1995); R. Baier, Yu. L. Dok- 
shitzer, A. Mueller, S. Peigne and D. Schiff, Nucl. Phys. 
B484, 265 (1997). 

B. G. Zakharov, JETP Letters, 63, 952 (1996); 65, 516 
(1997). 



M. Gyulassy, P. Levai, I. Vitev, iep-ph/9907461 



X.-N. Wang, Z. Huang and I. Sarcevic, Phys. Rev. Lett. 
77, 231 (1996); X.-N. Wang and Z. Huang, Phys. Rev. G 
55, 3047 (1997). 

X.-N. Wang, Phys. Rev. C 58, 2321 (1998); Phys. Rev. 
C 61, 064910 (2000). 

M. Luo, J.-W. Qiu and G. Sterman, Phys. Lett. B279, 
377 (1992); Phys. Rev. D49, 4493 (1994); Phys. Rev. 
D50, 1951 (1994). 

R. Field, Application of Perturbative QCD (Addison- 
Wesley, USA, 1989). 

X.-F. Guo, Phys. Rev. D 58, 114033(1998); X.-F. Guo 

and J. Qiu, Phys. Rev. D 61, 096003 (2000). 

J.-W. Qiu and G. Sterman, Nucl. Phys. B353 (1991) 105, 

137. 

X.-N. Wang and X. Guo, to be published. 
G. Aharelh and G. Parisi, Nucl. Phys. B126, 298 (1977). 
A. Mueller and J. Qiu, Nucl.Phys. B268, 427 (1986). 
J. Binnewies, B. A. Kniehl and G. Kramer, Z. Phys. C65, 
471 (1995). 

X.-F. Guo, Phys. Rev. D 58, 36001(1998). 
X. Guo, J. Qiu and X. Zhang, |hep-ph/991147e . 
X.-N. Wang, Phys. Lett. B485,157 (2000). 



4 



